Deaths with COVID-19 and from all-causes following first-ever SARS-CoV-2 infection in individuals with preexisting mental disorders: A national cohort study from Czechia

Background Evidence suggests reduced survival rates following Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection in people with preexisting mental disorders, especially psychotic disorders, before the broad introduction of vaccines. It remains unknown whether this elevated mortality risk persisted at later phases of the pandemic and when accounting for the confounding effect of vaccination uptake and clinically recorded physical comorbidities. Methods and findings We used data from Czech national health registers to identify first-ever serologically confirmed SARS-CoV-2 infections in 5 epochs related to different phases of the pandemic: 1st March 2020 to 30th September 2020, 1st October 2020 to 26th December 2020, 27th December 2020 to 31st March 2021, 1st April 2021 to 31st October 2021, and 1st November 2021 to 29th February 2022. In these people, we ascertained cases of mental disorders using 2 approaches: (1) per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes for substance use, psychotic, affective, and anxiety disorders; and (2) per ICD-10 diagnostic codes for the above mental disorders coupled with a prescription for anxiolytics/hypnotics/sedatives, antidepressants, antipsychotics, or stimulants per the Anatomical Therapeutic Chemical (ATC) classification codes. We matched individuals with preexisting mental disorders with counterparts who had no recorded mental disorders on age, sex, month and year of infection, vaccination status, and the Charlson Comorbidity Index (CCI). We assessed deaths with Coronavirus Disease 2019 (COVID-19) and from all-causes in the time period of 28 and 60 days following the infection using stratified Cox proportional hazards models, adjusting for matching variables and additional confounders. The number of individuals in matched-cohorts ranged from 1,328 in epoch 1 to 854,079 in epoch 5. The proportion of females ranged from 34.98% in people diagnosed with substance use disorders in epoch 3 to 71.16% in individuals diagnosed and treated with anxiety disorders in epoch 5. The mean age ranged from 40.97 years (standard deviation [SD] = 15.69 years) in individuals diagnosed with substance use disorders in epoch 5 to 56.04 years (SD = 18.37 years) in people diagnosed with psychotic disorders in epoch 2. People diagnosed with or diagnosed and treated for psychotic disorders had a consistently elevated risk of dying with COVID-19 in epochs 2, 3, 4, and 5, with adjusted hazard ratios (aHRs) ranging from 1.46 [95% confidence intervals (CIs), 1.18, 1.79] to 1.93 [95% CIs, 1.12, 3.32]. This patient group demonstrated also a consistently elevated risk of all-cause mortality in epochs 2, 3, 4, and 5 (aHR from 1.43 [95% CIs, 1.23, 1.66] to 1.99 [95% CIs, 1.25, 3.16]). The models could not be reliably fit for psychotic disorders in epoch 1. People diagnosed with substance use disorders had an increased risk of all-cause mortality 28 days postinfection in epoch 3, 4, and 5 (aHR from 1.30 [95% CIs, 1.14, 1.47] to 1.59 [95% CIs, 1.19, 2.12]) and 60 days postinfection in epoch 2, 3, 4, and 5 (aHR from 1.22 [95% CIs, 1.08, 1.38] to 1.52 [95% CIs, 1.16, 1.98]). Cases ascertained based on diagnosis of substance use disorders and treatment had increased risk of all-cause mortality in epoch 2, 3, 4, and 5 (aHR from 1.22 [95% CIs, 1.03, 1.43] to 1.91 [95% CIs, 1.25, 2.91]). The models could not be reliably fit for substance use disorders in epoch 1. In contrast to these, people diagnosed with anxiety disorders had a decreased risk of death with COVID-19 in epoch 2, 3, and 5 (aHR from 0.78 [95% CIs, 0.69, 0.88] to 0.89 [95% CIs, 0.81, 0.98]) and all-cause mortality in epoch 2, 3, 4, and 5 (aHR from 0.83 [95% CIs, 0.77, 0.90] to 0.88 [95% CIs, 0.83, 0.93]). People diagnosed and treated for affective disorders had a decreased risk of both death with COVID-19 and from all-causes in epoch 3 (aHR from 0.87 [95% CIs, 0.79, 0.96] to 0.90 [95% CIs, 0.83, 0.99]), but demonstrated broadly null effects in other epochs. Given the unavailability of data on a number of potentially influential confounders, particularly body mass index, tobacco smoking status, and socioeconomic status, part of the detected associations might be due to residual confounding. Conclusions People with preexisting psychotic, and, less robustly, substance use disorders demonstrated a persistently elevated risk of death following SARS-CoV-2 infection throughout the pandemic. While it cannot be ruled out that part of the detected associations is due to residual confounding, this excess mortality cannot be fully explained by lower vaccination uptake and more clinically recorded physical comorbidities in these patient groups.

Anatomical Therapeutic Chemical (ATC) classification codes.We matched individuals with preexisting mental disorders with counterparts who had no recorded mental disorders on age, sex, month and year of infection, vaccination status, and the Charlson Comorbidity Index (CCI).We assessed deaths with Coronavirus Disease 2019 (COVID- 19) and from allcauses in the time period of 28 and 60 days following the infection using stratified Cox proportional hazards models, adjusting for matching variables and additional confounders.The number of individuals in matched-cohorts ranged from 1,328 in epoch 1 to 854,079 in epoch 5.The proportion of females ranged from 34.98% in people diagnosed with substance use disorders in epoch 3 to 71.16% in individuals diagnosed and treated with anxiety disorders in epoch 5.The mean age ranged from 40.97 years (standard deviation [SD] = 15.69 years) in individuals diagnosed with substance use disorders in epoch 5 to 56.04 years (SD = 18.37 years) in people diagnosed with psychotic disorders in epoch 2. People diagnosed with or diagnosed and treated for psychotic disorders had a consistently elevated risk of dying with COVID-19 in epochs 2, 3, 4, and 5, with adjusted hazard ratios (aHRs) ranging from 1.46 [95% confidence intervals (CIs), 1.18, 1.79] to 1.93 [95% CIs, 1.12, 3.32].This patient group demonstrated also a consistently elevated risk of all-cause mortality in epochs 2, 3, 4, and 5 (aHR from 1.43 [95% CIs, 1.23, 1.66] to 1.99 [95% CIs, 1. 25, 3.16]).The models could not be reliably fit for psychotic disorders in epoch 1.People diagnosed with substance use disorders had an increased risk of all-cause mortality 28 days postinfection in epoch 3, 4, and 5 (aHR from 1.30 [95% CIs, 1.14, 1.47] to 1.59 [95% CIs, 1.19, 2.12]) and 60 days postinfection in epoch 2, 3, 4, and 5 (aHR from 1.22 [95% CIs, 1.08, 1.38] to 1.52 [95% CIs, 1.16, 1.98]).Cases ascertained based on diagnosis of substance use disorders and treatment had increased risk of all-cause mortality in epoch 2, 3, 4, and 5 (aHR from 1.22 [95% CIs, 1.03, 1.43] to 1.91 [95% CIs, 1.25, 2.91]).The models could not be reliably fit for substance use disorders in epoch 1.In contrast to these, people diagnosed with anxiety disorders had a decreased risk of death with COVID-19 in epoch 2, 3, and 5 (aHR from 0.78 [95% CIs, 0.69, 0.88] to 0.89 [95% CIs, 0.81, 0.98]) and all-cause mortality in epoch 2, 3, 4, and 5 (aHR from 0.83 [95% CIs, 0.77, 0.90] to 0.88 [95% CIs, 0.83, 0.93]).People diagnosed and treated for affective disorders had a decreased risk of both death with COVID-19 and from all-causes in epoch 3 (aHR from 0.87 [95% CIs, 0.79, 0.96] to 0.90 [95% CIs, 0.83, 0.99]), but demonstrated broadly null effects in other epochs.Given the unavailability of data on a number of potentially influential confounders, particularly body mass index, tobacco smoking status, and socioeconomic status, part of the detected associations might be due to residual confounding.

Conclusions
People with preexisting psychotic, and, less robustly, substance use disorders demonstrated a persistently elevated risk of death following SARS-CoV-2 infection throughout the pandemic.While it cannot be ruled out that part of the detected associations is due to residual confounding, this excess mortality cannot be fully explained by lower vaccination uptake and more clinically recorded physical comorbidities in these patient groups.
publicly deposited.For sharing the data, access must be granted by IHIS, the curator of Czech national health registers.Access to data can be requested through the procedure described at the web page of IHIS: https://www.uzis.cz/index-en.php?pg=contact-provision-of-information.The full analytical code is available at a dedicated GitHub repository: https://github.com/libpot/COVID_Mortality_Preexisting_Mental_Disorders.HM and JJ had full access to all registers and take responsibility for the integrity of the data export.TF, LP, HM, KM, and JJ had full access to all exported data.TF and LP responsibility for the accuracy of the data analysis.

Why was this study done?
• Existing research has demonstrated consistently elevated risk of death with Coronavirus Disease 2019 (COVID- 19) or all-cause mortality in people with preexisting psychotic and substance use disorders following a Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection.
• The evidence on people with preexisting affective and anxiety disorders is broadly consistent with increased mortality risk; however, multiple studies demonstrated null effects.
• To the best of our knowledge, no study has used national data covering almost all inpatient and outpatient settings, including primary care, and laboratory-confirmed SARS-CoV-2 infections to investigate whether this elevated mortality risk was present throughout the pandemic, including its later phases, and when robustly accounting for the confounding effect of vaccination uptake and clinically recorded physical comorbidities.
What did the researchers do and find?
• Using Czech national, whole population, all healthcare encompassing register-based data, individuals with preexisting psychotic and, less consistently, substance use disorders had increased risk of death with COVID-19 and all-cause mortality, including at the later phases of the pandemic.
• People with preexisting anxiety disorders had decreased risk of death with COVID-19 and all-cause mortality in multiple epochs, whereas people with preexisting affective disorders demonstrated broadly null effects throughout the pandemic.
• These associations could not be fully explained by differences in vaccination uptake or clinically recorded physical comorbidities.

What do these findings mean?
• The consistently lower survival in people with preexisting psychotic and substance use disorders aligns with existing evidence on fatal health inequalities in these patient groups.
• Systemic efforts are needed to fully reverse the risk attributable to long-term, structural processes affecting health of people with psychotic and substance use disorders.
• The main limitation of the present study was its inability to fully control for a number of characteristics, particularly body mass index, tobacco smoking status, and socioeconomic status, that might confound the associations between mental disorders and mortality: future studies should explore these associations while accounting for these confounders.

Introduction
Evidence before the outbreak of the Coronavirus Disease 2019 (COVID-19) pandemic showed that people with mental disorders have a higher risk of developing a wide range of physical health conditions relative to their counterparts without these disorders [1][2][3] as well as higher mortality rates and shorter life expectancies than the general population [2,[4][5][6][7].Worse general health, often associated with lower socioeconomic status and lifestyle risk factors (e.g., smoking) could contribute to an increased risk of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection and potentially lower survival following the infection in these people.
While this evidence suggests lower survival following SARS-CoV-2 infection in people with preexisting mental disorders, to the best of our knowledge, no study has used national data covering almost all inpatient and outpatient settings, including primary care, and laboratoryconfirmed SARS-CoV-2 infections to investigate whether this elevated mortality risk persisted at later phases of the pandemic and when robustly accounting for the confounding effect of vaccination uptake and clinically recorded physical comorbidities.
In the present study, we used national, whole population, all healthcare encompassing register-based data to investigate the risk of death with COVID-19 and from all-causes following first-ever laboratory-confirmed infection with SARS-CoV-2 in individuals with preexisting mental disorders compared with matched counterparts without mental disorders at 5 distinct pandemic phases.By performing matching on vaccination status and clinically recorded physical comorbidities, we aimed to explore associations not confounded by differences between people with and without preexisting mental disorders on these characteristics.

Methods
The research questions and the analytical plan were preregistered at Open Science Framework before data analyses started [24]: any deviations from the plan are described in S1 Methods.This study was reported as per the Reporting of studies Conducted using Observational Routinely collected health Data (RECORD) Statement (see S1 RECORD Checklist).

Setting
Mental health care in Central and Eastern European region relies on large psychiatric hospitals [25][26][27].Considering Czechia in particular, more than 50% of its mental health budget is allocated to inpatient services [28], with the majority of inpatient care provided in outdated psychiatric hospitals [29].However, Czechia has launched its mental health reform in 2013, with its initial main goals focusing on deinstitutionalization.This entails the expansion of community-based services, alongside a reduction in long-term inpatient beds and complemented by educational, destigmatization and other implementation programs aimed at improving the quality of care and overall quality of life of people with psychiatric conditions [27,30].
The first wave of the pandemic in Czechia lasted roughly from 1st March 2020 to 30th September 2020, with a State of Emergency being in place from 12th March 2020 to 17th May 2020.The first wave resulted in 70,968 incident infections [31].The second wave of the pandemic lasted approximately from 1st October 2020 to 31th March 2020, with a State of Emergency imposed from 5th October 2020 to 11th April 2021.The second wave led to 1,482,727 incident infections [31].The Czech National Vaccination Strategy was launched in December 2020, with preexisting mental disorders not considered as reason for priority inoculation [32].The period of post-second wave lasted roughly from 1st April 2021 to 31st October 2021, resulting in 227,827 incident infections [31].Then, the 2021 to 2022 winter wave lasted from approximately 1st November 2021 to 28th February 2022, with a State of Emergency imposed from 26th November to 25th December, and led to 1,820,446 incident infections [31].

Ethics statement
This study was approved by the Ethics Committee of the National Institute of Mental Health (approval number 176/21).

Data
We used data from the National Registry of Reimbursed Health Services (NRRHS), part of the National Health Information System (NHIS), covering inpatient and outpatient services, including primary care, as well as prescription medications.The register covers nearly the entire Czech population (approximately 10.7 million inhabitants).The records are created by health professionals who complete information on diagnosis (primary and secondary diagnoses) as per the International Classification of Diseases 10th Revision (ICD-10), date (for inpatient settings, admission, and discharge date), Anatomical Therapeutic Chemical (ATC) classification codes for prescription medications (with the exception of common medications administered in inpatient settings), and basic sociodemographic information such as age, sex, and region of permanent residency.Additionally, we used data from the Information System of Infectious Diseases (ISID) covering nationwide testing for SARS-CoV-2 and COVID-19 vaccination status.Furthermore, we used data from the register of all-cause mortality, containing information on the date of death, the ICD-10 cause(s), and, if applicable, the external cause(s) of death.All 3 registers can be interlinked using a common unique identifier and are maintained by the state-funded Institute of Health Information and Statistics of Czechia (IHIS).Data in the NHIS are collected in accordance with Act No. 372/2011 Coll., on health services and conditions of their provision, while ISID data are collected in accordance with Act No. 258/2000 Coll., on public health protection.Due to this legal mandate, the retrospective analyses of data in these registries did not require informed consents from participants.
We retrieved all individuals aged 10 or above-the earliest plausible onset age of the studied mental disorders [33]-with first-ever laboratory-confirmed SARS-CoV-2 infection occurring in 5 epochs: 1. 1st March 2020 to 30th September 2020, the first wave of the pandemic.

Exposure
In individuals with first-ever laboratory-confirmed SARS-CoV-2 infection, we used 2 approaches to ascertain cases.The first approach relied on identifying the occurrence of diagnosis per ICD-10 codes for (1) substance use disorders (F1); (2) psychotic disorders (F2); (3) affective disorders (F3); and (4) anxiety disorders (F4) in the period of 5 years prior to the date of infection (see details in S1 Methods).We considered the occurrence of at least one of the above codes as any mental disorder.We established the occurrence of each of the mental disorders separately.We considered an individual to have a diagnosis when the given ICD-10 code was listed on a record in either inpatient (primary diagnosis, considered from discharge date) or any outpatient setting.Conversely, the unexposed cohort included individuals who had no such occurrence in the period of 5 years before the date of their SARS-CoV-2 infection.
The second approach entailed establishing whether an individual was prescribed psychopharmaceuticals at least once in the period of 5 years prior to the date of SARS-CoV-2 infection, in addition to occurrence of diagnosis (inpatient or any outpatient setting) for a given ICD-10 code.We considered the prescription per ATC codes of any anxiolytics/hypnotics/sedatives (N05B, N05C), antidepressants (N06A), antipsychotics (N05A), or stimulants (N06B).Conversely, the unexposed cohort included individuals who had no diagnosis of a mental disorder and no prescription of any psychopharmaceutical in the period of 5 years before the date of their infection.
We used the 2 ascertainment approaches to investigate the consistency of estimates across different exposure definitions: broadly consistent results between these would increase the confidence in the robustness of inferences.

Control of confounding
In our identification and selection of potential confounders, we followed the "disjunctive cause criterion," in which one controls for covariates that are causes of the exposure or causes of the outcome or causes of both [34,35].

Matching
In the first 2 epochs, we matched on age, sex, month and year of infection as well as the Charlson Comorbidity Index (CCI) [36].In the 3 subsequent epochs, we matched on age, sex, month and year of infection, vaccination status, and the CCI.Since vaccination does not confer an immediate protection, we did not consider vaccinations that were administered 14 or less days before the infection.For example, when an individual received the first dose of a twodose regimen more than 14 days before the infection, and the second dose 14 or less days before the infection, we considered them as having received the first dose at the time of the infection.The CCI referred to the period of 5 years before the date of the SARS-CoV-2 infection and was coded as 0, 1, 2, 3, and 4 or more comorbidities.Each exposed individual was matched with up to 5 unique unexposed counterparts.Some people with preexisting mental disorders had no matching counterparts (in each cohort <10%, see details in S1-S5 Tables); we excluded these unmatched individuals from the respective analyses.

Additional confounders
To further reduce the level of unaccounted for confounding, we adjusted for region of permanent residency, overall number of contacts with inpatient services and overall number of contacts with outpatient services (disregarding contacts related to the exposure), and prescription medications (see the detailed list with ATC codes in S1 Methods).We considered prescription of each of the medications or treatment administration in the period of 1 year prior to the SARS-CoV-2 infection separately.The number of contacts with the healthcare system referred to the period of 5 years before the date of the SARS-CoV-2 infection.For details, see the proposed directed acyclic graph in S1 Fig in S1 File.

Outcome
We considered (1) deaths with COVID-19 (ICD-10 codes U071 and U072 listed as a cause of death on the death certificate); and (2) all-cause mortality occurring in the period of (1) 28 days; and (2) 60 days after a positive test for SARS-CoV-2.These cutoffs are based on Public Health England's analysis that showed that 88% and 96% of deaths occurred within 28 and 60 days of a positive test, respectively [37].

Statistical analysis
Following descriptive analysis, we used stratified Cox proportional hazards models to assess the risks of deaths with COVID-19 and from all-causes in individuals with preexisting mental disorders compared with matched counterparts without such disorders, separately for each studied mental disorder and epoch.Each stratum consisted of 1 person with preexisting mental disorder and up to 5 matched counterparts.Time-to-event was expressed in days.In models investigating the risk of death with COVID-19, we considered death due to any other cause as competing risk, and the affected individuals were censored.We fitted models adjusting for confounders, with the CCI used as a continuous measure.The results were expressed as hazard ratios (HRs) with 95% confidence intervals (95% CIs).We tested the proportionality assumption using Schoenfeld residuals; in some instances, the assumption was violated, we therefore interpreted the HRs as weighted averages of the time-varying HRs over the entire follow-up period [38].In line with the statement from the American Statistical Association on p-values [39], we present effect sizes with 95% CIs throughout the manuscript.However, we provide pvalues as complementary information in Supplementary Results.All analyses were conducted in R statistical programming language (version 4.2.2) [40], using the libraries survival (version 3.5-5) and EValue (version 4.1.3)[41].

Sensitivity analyses
Having a history of a mental disorder might influence the risk of being tested for SARS-CoV-2 infection; thus, restricting the analysis to individuals who had a positive test might lead to collider bias [42].To examine potential presence of collider bias, we conducted negative control exposure analyses by assessing the associations between the characteristics that are expected to be unrelated to the outcome and the outcome itself within the chosen cohorts [42].To do so, we considered the occurrence of ( 1) migraine (ICD-10 code G43); (2) fracture of forearm (ICD-10 code S52); (3) acne (ICD-10 code L70); (4) mild allergies (ICD-10 codes J301, L500, and L23); and (5) transport accidents (ICD-10 codes V01-V99) in the time period of 5 years prior to the positive SARS-CoV-2 test.We assessed the occurrence of each of these separately.Then, we fitted stratified Cox proportional hazards models with the above characteristics being the exposures, while using the confounders and outcomes from the main analysis.We reported the total number and proportion of non-null tests, with the theoretical maximum being 40 (or 100%) per 1 epoch-mental disorder combination (i.e., 2 case ascertainment definitions X 5 negative control exposures X 4 outcomes).Since some of our cohorts were considerably large (i.e., anxiety disorders in epochs 2 to 5), it would be possible to have non-null results even if the effect sizes were negligible [43].Thus, we complementarily provided averaged HRs across the 40 tests per 1 epoch-mental disorder combination.Proportion of non-null tests closer to 0% and/or averaged HRs closer to 1 would suggest the absence of collider bias.
To assess the level of potential unmeasured confounding, we calculated E-values for each of our regression model where the results were not consistent with a null effect.The E-value indicate the strength of association-here expressed in hazard ratio-an unmeasured confounder, or set of confounders, would need to have with both the exposure and the outcome to nullify the association between the exposure and the outcome observed in the model [44].E-values closer to 1 indicate lower confidence in the results not being due to residual confounding [44].

Results
The number of individuals in matched-cohorts ranged from  1 and 2 and those with additional confounders in S6 and S7 Tables.

Risk of death with COVID-19 in people with preexisting mental disorders
In the models adjusting for all considered confounders, including vaccination uptake and clinically recorded physical comorbidities, people diagnosed with or diagnosed and treated for psychotic disorders had an elevated risk of death with COVID-19 in epochs 2, 3, 4, and 5, both 28 and 60 following SARS-CoV-2 infection.The models could not be reliably fit for psychotic disorders in epoch 1.Those diagnosed with substance use disorders had an increased risk of death with COVID-19 28 days postinfection in epoch 3 and 4 and 60 days postinfection in epoch 3. Cases ascertained based on diagnosis of substance use disorders and treatment by psychopharmaceuticals had an elevated risk of death with COVID-19 in epoch 3, both 28 and 60 days following infection.The models could not be reliably fit for substance use disorders in epoch 1, and the remaining ones were consistent with a null effect.
In contrast, people diagnosed with or diagnosed and treated for anxiety disorders had a decreased risk of death with COVID-19 in epoch 2, 3, and 5, both 28 and 60 days postinfection.The remaining models for anxiety disorders were consistent with a null effect.Additionally, people diagnosed and treated for affective disorders had a decreased risk of death with COVID-19 in epoch 3, both 28 and 60 days postinfection, but all other models involving affective disorders were broadly consistent with a null effect.The results for any studied mental disorder were-regardless of case ascertainment definition-consistent with a null effect in all epochs.For detailed results, see Figs 1 and S2-101 in S1 File and S8-S10 Tables.

Risk of all-cause mortality in people with preexisting mental disorders
In the models adjusting for all considered confounders, including vaccination uptake and clinically recorded physical comorbidities, people diagnosed with or diagnosed and treated for psychotic disorders were more likely to die in epochs 2, 3, 4, and 5, both 28 and 60 days postinfection.The models could not be reliably fit for psychotic disorders in epoch 1.In those diagnosed with substance use disorders, there was an elevated risk of all-cause mortality 28 days postinfection in epoch 3, 4, and 5 and 60 days postinfection in epoch 2, 3, 4, and 5. Cases ascertained based on diagnosis of substance use disorders and treatment had increased risk of all-cause mortality in epoch 2, 3, 4, and 5, both 28 and 60 days postinfection.The models could not be reliably fit for substance use disorders in epoch 1, and the remaining ones were consistent with a null effect.

Sensitivity analyses
In negative control exposure analyses, the proportion of non-null tests did not exceed 15% for substance use disorders, 15% for psychotic disorders, 20% for affective disorders, 40% for For models not consistent with a null-effect, the E-values ranged from 1.71 to 3.22 for substance use disorders, from 2.21 to 3.39 for psychotic disorders, from 1.50 to 1.88 for anxiety disorders, and from 1.45 to 1.55 for affective disorders.E-values further away from 1 increase the confidence that the detected associations are not due to unaccounted for confounding.See details in S15 and S16 Tables.

Discussion
Using Czech national health register data, we demonstrated that people with preexisting psychotic disorders were more likely to die with COVID-19 or due to any cause following SARS-CoV-2 infection throughout the pandemic.We demonstrated less robust associations for deaths with COVID-19 in people with substance use disorders but they had a consistent and sustained elevated risk of all-cause mortality following SARS-CoV-2 infection.The 2 exposure definitions produced broadly consistent results across each epoch-mental disorder combination.This detected excess mortality is not fully explicable by differences in vaccination uptake or clinically recorded physical comorbidities between people with and without preexisting substance use and psychotic disorders.Separately, people with anxiety disorders demonstrated decreased risk of death with COVID-19 and from all-causes in multiple epochs, whereas the risk in people with affective disorders was broadly consistent with a null effect throughout the pandemic.
Our findings are in line with existing evidence showing elevated mortality risk in people with psychotic [10][11][12][13][14][15][16][17] and substance use disorders [18][19][20][21]; however, we demonstrated that these health disparities were consistent throughout the pandemic and persisted at its later phases too.Robust control for vaccination uptake or clinically recorded somatic comorbidity in our study did not reverse the increased risk of death.
These results, broadly consistent with a citywide study from the United Kingdom [45], demonstrate vulnerability in these patient groups that cannot be fully explained by differences in vaccination uptake or clinically recorded physical comorbidities.This suggest that other individual and structural factors might be responsible for the detected outcomes.Inequalities in access to healthcare and differences in the quality of care received cannot be discounted as contributing to the excess mortality.Previous studies showed that these patient groups may face delayed diagnosis [46,47] if physical health conditions are recognized at all [48]; such suboptimal episodes of care may be related to, among other things, incorrectly attributing the symptoms of somatic conditions to mental disorders [49].Thus, people with psychotic and substance use disorders potentially have more severe and insufficiently addressed or clinically unrecognized physical comorbidities that contributed to lower survival.
People with lower socioeconomic status were less likely to be tested for SARS-CoV-2 infection [50][51][52], but more likely to experience delayed test results [53].Both substance use and psychotic disorders are negatively associated with socioeconomic status [54][55][56][57].It is plausible that SARS-CoV-2 infection in these people was recognized comparatively late and that this potentially adversely influenced the therapeutic response.
Further, negative health behaviors such as smoking tobacco and suboptimal nutrition and physical activity are common in these patient groups [58,59] and may have contributed to worse prognosis.Pharmacological treatments for psychotic disorders are known to contribute to metabolic disturbances [60], interact with or limit the use of treatment for somatic conditions, thus potentially contributing to lower survival also post-SARS-CoV-2 infection.However, pre-pandemic research has demonstrated that the use of antipsychotics is associated with decreased risk of all-cause mortality in people with schizophrenia [61], with no differences between concomitant use of several ones compared with monotherapy [62].Alternatively, lower adherence to prescription medications [63] that would contribute to worsened overall health at the baseline, cannot be ruled out.
Overall, the existing evidence and our results suggest the presence of fatal but largely tractable health inequalities in these patient groups.Wholesale system approaches, predicated on evidence-based policy changes and, ideally, combined with evaluation, are required to address the multilayered factors behind these fatal health inequalities.In a future pandemic or other health emergency, substance use and psychotic disorders need to be considered as a specific vulnerability factor beyond liability to the health threat, itself.
Regarding other psychiatric conditions, our study further showed decreased mortality risk in people with anxiety disorders in multiple epochs, contrasting findings of existing studies that demonstrated elevated risks [11,13,16,22] or null effects [15,18].Multiple factors might be responsible for these differences, including the scope of data, definition of cases and the comparison groups, as well as analytical approaches.In 2 French studies, for example, anxiety disorders were identified through a hospital register [18,22], while our study included cases in all healthcare settings, likely capturing less severe cases with better outcomes.In a UK Biobank study, narrower definitions of anxiety-and stress-related disorders were used [13], focusing on more severe conditions.The study also did not use a matched-cohort design [13], making it difficult to rule out that adjusting for clinical, sociodemographic, and behavioral confounders in regression models still led to residual confounding due to covariate imbalance.Similarly, in a QResearch database study, the authors did not implement a matched-cohort design and also did not condition on positive test for a SARS-CoV-2 infection, hence, comparisons involved different populations [16].
Evidence shows that anxiety symptoms may have been linked to higher endorsement of preventive measures by heightened contamination fear [64,65].Together with our findings of decreased mortality risk, this may suggest that excessive preoccupation with COVID-19related events may have facilitated early detection and improved the health outcomes in some people with anxiety disorders.
The broadly null effect we detected for mortality in people with affective disorders was reported before also elsewhere [15]; however, other studies demonstrated an elevated risk [11,14,16,23].As with anxiety disorders, the factors responsible for these diverging results are most likely multiple and include breadth of data (all healthcare settings versus specific settings), case definitions (all affective disorders diagnostic codes versus specific diagnoses of affective disorders), and analytical choices (matched-cohorts versus unmatched-cohorts, conditioning on positive SARS-CoV-2 infection tests versus not conditioning on these).
Strengths included the use of national, whole population, fully standardized data on SARS-CoV-2 infection status, healthcare utilization, and mortality.Next, we investigated the robustness of our results through using multiple definitions of the exposure and the outcome.Then, we used negative control exposures to explore the potential presence of collider bias and E-values to establish what level of unaccounted for confounding would explain away the observed associations.
This study has also some limitations.First, we used broad diagnostic categories of mental disorders; thus disregarding the diagnostic heterogeneity in these, including the differing levels of severity (e.g., bipolar disorder in affective disorders).Second, we matched on key sociodemographic and clinical covariates and subsequently adjusted for a wide range of additional health-related confounders; however, we did not control for a number of previously identified influential clinical, sociodemographic, and behavioral confounders [66].In particular, we had no information on body mass index, smoking status, and socioeconomic status per se; however, these are likely to influence an individual's overall health [67][68][69][70][71][72], which we considered by controlling for a comorbidity index, overall number of contacts with inpatient and outpatient services, and prescription medications, including antihypertensives and statins that would be administered for conditions commonly present in people who smoke or who are obese [73][74][75][76].These steps likely reduced the level of confounding due to these covariates; however, we cannot rule out that part of the detected associations is still due to residual confounding, both related to these known and directly unmeasured and potential unknown confounders, with the most plausible direction of bias being the overestimation of true effects.Third, we were able to match the vast majority of people with mental disorders with counterparts without mental disorders; however, "bias due to incomplete matching" cannot be ruled out [77].In particular, the unmatched individuals with mental disorders tended to be, on average, older and have a higher number of comorbidities.Since these individuals can be expected to have the worse outcomes postinfection, the most plausible direction of bias seems to be the underestimation of true effects.Fourth, while clinically recorded physical comorbidities are among the factors most strongly associated with worse prognosis following a SARS-CoV-2 infection [78], we cannot rule out that in some individuals with mental disorders, they would act as a mediator instead of a confounder, thus raising the possibility of overadjustment bias [79].Fifth, both the curator of the data, IHIS, and the insurance companies who use them to reimburse service providers employ mechanisms to ensure the validity of data; however, all diagnoses used in this study have not been fully validated yet.Thus, under-registration and/or errors in diagnoses coding cannot be ruled out.Sixth, some of our analyses included considerably few individuals, leading to profound uncertainty in our estimates.Seventh, we did not investigate survival following SARS-CoV-2 re-infections in people with preexisting disorders.Eighth, we did not investigate the responses to medications and/or other treatment modalities following the infection itself in people with preexisting mental disorders.Ninth, we did not consider the outcomes of people with multiple psychiatric conditions.Last, the follow-up period following infection was short; however, we had no information on emigration status, and we cannot rule out that some individuals were lost to follow-up.
People with preexisting psychotic, and, less robustly, substance use disorders demonstrated persistently lower survival following SARS-CoV-2 infection throughout the pandemic.While we cannot rule out that part of the detected associations is due to residual confounding, the consistently increased vulnerability beyond vaccination uptake and clinically recorded physical comorbidity aligns with existing evidence on fatal health inequalities in these patient groups and underlines the importance of implementing systemic efforts to fully reverse these.To at least reduce these disparities, it must be assured that these patient groups are included in future vaccination campaigns.

2 .
1st October 2020 to 26th December 2020, the second wave of the pandemic before the initiation of the national vaccination program.

3 .
27th December 2020 to 31st March 2021, the beginning of the national vaccination program to the end of the second wave of the pandemic.4.1st April 2021 to 31st October 2021, the post-second wave period.5. 1st November 2021 to 28th February 2022, the 2021 to 2022 winter wave.

10 )
diagnostic codes.S4 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with any mental disorder in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S5 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with any mental disorder in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S6 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with any mental disorder in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S7 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with any mental disorder in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S8 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with any mental disorder in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S9 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with any mental disorder in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S10 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with any mental disorder in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S11 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with any mental disorder in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S12 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with substance use disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S13 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with substance use disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S14 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with substance use disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S15 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with substance use disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S16 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with substance use disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S17 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with substance use disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S18 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with substance use disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S19 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with substance use disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S20 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with substance use disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S21 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with substance use disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S22 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with psychotic disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S23 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with psychotic disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S24 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with psychotic disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S25 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with psychotic disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S26 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with psychotic disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S27 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with psychotic disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S28 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with psychotic disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S29 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with psychotic disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S30 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with psychotic disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S31 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with psychotic disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S32 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with affective disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S33 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with affective disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S34 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with affective disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S35 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with affective disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S36 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with affective disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S37 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with affective disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S38 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with affective disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S39 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with affective disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S40 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with affective disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S41 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with affective disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S42 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with anxiety disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S43 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with anxiety disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S44 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with anxiety disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S45 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with anxiety disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S46 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with anxiety disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S47 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with anxiety disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S48 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with anxiety disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S49 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with anxiety disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S50 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with anxiety disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S51 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 28 days in people with anxiety disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S52 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with any mental disorder in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S53 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with any mental disorder in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S54 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with any mental disorder in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S55 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with any mental disorder in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S56 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with any mental disorder in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S57 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with any mental disorder in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S58 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with any mental disorder in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S59 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with any mental disorder in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S60 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with any mental disorder in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S61 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with any mental disorder in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S62 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with substance use disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S63 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with substance use disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S64 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with substance use disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S65 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with substance use disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S66 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with substance use disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S67 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with substance use disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S68 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with substance use disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S69 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with substance use disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S70 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with substance use disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S71 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with substance use disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S72 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with psychotic disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S73 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with psychotic disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S74 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with psychotic disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S75 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with psychotic disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S76 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with psychotic disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S77 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with psychotic disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S78 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with psychotic disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S79 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with psychotic disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S80 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with psychotic disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S81 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with psychotic disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S82 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with affective disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S83 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with affective disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S84 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with affective disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S85 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with affective disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S86 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with affective disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S87 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with affective disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S88 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with affective disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S89 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with affective disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S90 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with affective disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S91 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with affective disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S92 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with anxiety disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S93 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with anxiety disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S94 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with anxiety disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S95 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with anxiety disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S96 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with anxiety disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S97 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with anxiety disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S98 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with anxiety disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S99 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with anxiety disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S100 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with anxiety disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S101 Fig. Cumulative probability with 95% confidence interval of death with COVID-19 up to 60 days in people with anxiety disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S102 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with any mental disorder in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S103 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with any mental disorder in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S104 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with any mental disorder in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S105 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with any mental disorder in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S106 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with any mental disorder in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S107 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with any mental disorder in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S108 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 28 days in people with any mental disorder in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S109 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with any mental disorder in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S110 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with any mental disorder in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S111 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with any mental disorder in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S112 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 28 days in people with substance use disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S113 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with substance use disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S114 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with substance use disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S115 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with substance use disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S116 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with substance use disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S117 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with substance use disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S118 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with substance use disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S119 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with substance use disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S120 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with substance use disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S121 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with substance use disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S122 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with psychotic disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S123 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with psychotic disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S124 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with psychotic disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S125 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with psychotic disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S126 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with psychotic disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S127 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with psychotic disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S128 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with psychotic disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S129 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with psychotic disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S130 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with psychotic disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S131 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 28 days in people with psychotic disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S132 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with affective disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S133 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 28 days in people with affective disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S134 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with affective disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S135 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with affective disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S136 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with affective disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S137 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with affective disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S138 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with affective disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S139 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with affective disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S140 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 28 days in people with affective disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S141 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with affective disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S142 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with anxiety disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S143 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with anxiety disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S144 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 28 days in people with anxiety disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S145 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with anxiety disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S146 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with anxiety disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S147 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with anxiety disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S148 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with anxiety disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S149 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 28 days in people with anxiety disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S150 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with anxiety disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S151 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 28 days in people with anxiety disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S152 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with any mental disorder in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S153 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with any mental disorder in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S154 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with any mental disorder in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S155 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with any mental disorder in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S156 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with any mental disorder in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S157 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with any mental disorder in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S158 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 60 days in people with any mental disorder in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S159 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with any mental disorder in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S160 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with any mental disorder in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S161 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with any mental disorder in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S162 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 60 days in people with substance use disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S163 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with substance use disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S164 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with substance use disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S165 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with substance use disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S166 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with substance use disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S167 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with substance use disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S168 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with substance use disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S169 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with substance use disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S170 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with substance use disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S171 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with substance use disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S172 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with psychotic disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S173 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with psychotic disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S174 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with psychotic disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S175 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with psychotic disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S176 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with psychotic disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S177 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with psychotic disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S178 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with psychotic disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S179 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with psychotic disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S180 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with psychotic disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S181 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 60 days in people with psychotic disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S182 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with affective disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S183 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 60 days in people with affective disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S184 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with affective disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S185 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with affective disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S186 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with affective disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S187 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with affective disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S188 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with affective disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S189 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with affective disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S190 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 60 days in people with affective disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S191 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with affective disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S192 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with anxiety disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S193 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with anxiety disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S194 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 60 days in people with anxiety disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S195 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with anxiety disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S196 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with anxiety disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes.S197 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with anxiety disorders in epoch 1, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S198 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with anxiety disorders in epoch 2, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S199 Fig. Cumulative probability with 95% confidence interval of allcause mortality up to 60 days in people with anxiety disorders in epoch 3, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S200 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with anxiety disorders in epoch 4, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.S201 Fig. Cumulative probability with 95% confidence interval of all-cause mortality up to 60 days in people with anxiety disorders in epoch 5, cases ascertained by diagnosis per the International Classification of Diseases 10th Revision (ICD-10) diagnostic codes coupled with prescription for psychopharmaceuticals per the Anatomical Therapeutic Chemical (ATC) classification codes.(ZIP)